Search for: All records

In this paper, we study the convergence analysis for a robust stochastic structure-preserving Lagrangian numerical scheme in computing effective diffusivity of time-dependent chaotic flows, which are modeled by stochastic differential equations (SDEs). Our numerical scheme is based on a splitting method to solve the corresponding SDEs in which the deterministic subproblem is discretized using a structure-preserving scheme while the random subproblem is discretized using the Euler-Maruyama scheme. We obtain a sharp and uniform-in-time convergence analysis for the proposed numerical scheme that allows us to accurately compute long-time solutions of the SDEs. As such, we can compute the effective diffusivity for time-dependent chaotic flows. Finally, we present numerical results to demonstrate the accuracy and efficiency of the proposed method in computing effective diffusivity for the time-dependent Arnold-Beltrami-Childress (ABC) flow and Kolmogorov flow in three-dimensional space. 

The G-equation is a well-known model for studying front propagation in turbulent combustion. In this paper, we develop an efficient model reduction method for computing regular solutions of viscous G-equations in incompressible steady and time-periodic cellular flows. Our method is based on the Galerkin proper orthogonal decomposition (POD) method. To facilitate the algorithm design and convergence analysis, we decompose the solution of the viscous G-equation into a mean-free part and a mean part, where their evolution equations can be derived accordingly. We construct the POD basis from the solution snapshots of the mean-free part. With the POD basis, we can efficiently solve the evolution equation for the mean-free part of the solution to the viscous G-equation. After we get the mean-free part of the solution, the mean of the solution can be recovered. We also provide rigorous convergence analysis for our method. Numerical results for viscous G-equations and curvature G-equations are presented to demonstrate the accuracy and efficiency of the proposed method. In addition, we study the turbulent flame speeds of the viscous G-equations in incompressible cellular flows. 

A mathematical model is proposed for shape evolution and locomotion of fish epidermal keratocytes on elastic substrates. The model is based on mechanosensing concepts: cells apply contractile forces onto the elastic substrate, while cell shape evolution depends locally on the substrate stress generated by themselves or external mechanical stimuli acting on the substrate. We use the level set method to study the behaviour of the model numerically, and predict a number of distinct phenomena observed in experiments, such as (i) symmetry breaking from the stationary centrosymmetric to the well-known steadily propagating crescent shape, (ii) asymmetric bipedal oscillations and travelling waves in the lamellipodium leading edge, (iii) response to remote mechanical stress externally applied to the substrate (tensotaxis) and (iv) changing direction of motion towards an interface with a rigid substrate (durotaxis).